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ABSTRACT 

We investigate the possibility of constructing models of i?-violating LQD Yukawa cou- 
plings using a single U(l) flavour-symmetry group and supermultiplet charge assignments 
that are compatible with the known hierarchies of quark and lepton masses. The mis- 
match of mass and current eigenstates inferred from the known charged-current mixing 
induces the propagation of strong phenomenological constraints on some R- violating cou- 
plings to many others. Applying these constraints, we look for flavour-symmetry models 
that are consistent with different squark-production hypotheses devised to explain the 
possible HERA large-Q 2 anomaly. The e + d — > t interpretation of the HERA data is 
accommodated relatively easily, at the price of postulating an extra parity. The e + s — > t 
interpretation of the events requires models to have only small (2,3) mixing in the down 
quark sector. The e + d — > c mechanism cannot be accommodated without large vio- 
lations of squark-mass universality, due to the very strong experimental constraints on 
i?-violating operators. We display a model in which baryon decay due to dangerous 
dimension-five operators is automatically suppressed. 
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1 Introduction 



Although the minimal supersymmetric Standard Model (MSSM) has dominated the phe- 
nomenological studies of supersymmetric signals it has long been known that the 
symmetries of the Standard Model allow additional dimension-four couplings which may 
lead to interesting baryon- and lepton- number- violating processes 0. These couplings 
are expected to be present in the low energy Lagrangian, unless forbidden by a symmetry 
such as R parity ||. The complete set of such terms in the superpotential is: 

XLiLjEk + X'LiQjDk + \"UiDjDk 

where the L(Q) are the left-handed lepton (quark) superfields, and the E,(D, U) are the 
corresponding right-handed fields. The symmetries of the model imply that there are 45 
operators in total. However, there are many experimental constraints on these operators 
and their combinations, of which the most stringent comes from proton stability and ex- 
cludes the simultaneous presence of certain products of LQD and U DD couplings || . In 
addition, experimental constraints from the non-observation of modifications to Standard 
Model processes, or of possible exotic processes, gives bounds for most of the operators 
and some combinations involving pairs of fermion generations. On the other hand, possi- 
ble strong limits on i?-parity violating interactions from cosmological arguments || can 
be avoided in various schemes pj, including the case of electroweak baryogenesis ||. 

The large number of i?-violating couplings complicates the systematic discussion of the 
phenomenological implications of these constraints. To date, most phenomenological 
analyses have assumed the dominance of a single operator, arguing that the Yukawa cou- 
plings of the Standard model display just such a property. In flavour-symmetry models, 
the dominant operator is naturally specified in the quark and lepton current basis. It 
is plausible to assume that mass mixing will induce non-zero coefficients for operators 
related to the dominant one. In Section 2 of this paper, we pursue this argument and com- 
pile the corresponding implications of some severe upper limits on particular i?-violating 
interactions. 

There have been many attempts to understand quark and lepton masses and the mixing 
angles between mass and current eigenstates using models for family symmetries [[J . Some 
of these reproduce successfully the qualitative features of fermion masses and mixings, 
and so provide plausible frameworks for analyzing the possible hierarchy of i?-violating 
interactions [|IIJ. In this paper we consider models based on a single U(l) family symme- 
try, with fermion charges constrained by the observed hierarchy of fermion masses and 
mixing angles [II]. Such models are discussed in Section 3, where problems arising from 



symmetric mass matrices and from constraints on products of operators are emphasized. 
As a specific application of this analysis, we look in Section 4 for models that might 



accommodate the proposed i?-violating interpretations |L2| of the possible HERA large- 
Q 2 anomaly 119, WM. Of the 45 operators mentioned earlier, 9 could in principle lead 
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to resonant squark production at HERA. Of these, only the X' 121 ,\' 131 and A' 132 cases 
survived an initial confrontation with other experimental constraints [12] The suggestion 
that the apparent HERA excess may be due to single sparticle production via some R- 



violating couplings may not be gaining support [15]. Nevertheless, our analysis gives 
an indication which of the proposed mechanisms may be compatible with U(l) family- 
symmetry models. Within the framework of the most symmetric schemes describing 
fermion masses, we find the bounds on R-violating couplings are so strong that such 
schemes do not lead to a significant excess of HERA events over the Standard Model 
prediction. However, in more general schemes we find that the A' 131 interpretation is 
easy to accommodate, whereas the A' 121 interpretation has difficulties with squark mass 
universality. We indicate how to construct a model consistent with the A' 132 interpretation, 
though we do not present a specific example. We also show how the structure of the quark 
and lepton mass matrices would be strongly constrained by the confirmation of such an 
.R-violating signal. 



2 R Violation and Family Symmetries 

In order to obtain a realistic form for the quark and lepton masses and mixing angles, 
it is necessary to have non-diagonal forms for the mass matrices in the current basis. 
Diagonalising the mass matrix then implies that the mass eigenstates are mixtures of the 
current eigenstates. Attempts to make sense of the pattern of fermion masses and mixing 
angles often start with a family symmetry in the current basis which, when exact, allows 
only the third generation of quarks and leptons to acquire mass. Spontaneous breaking 
of this symmetry then allows other entries of the mass matrix to be non-zero. If the 
breaking is weak, these entries will be small, offering an explanation for the observed 
hierarchy of fermion masses and mixing angles. 

If R parity is violated, such a symmetry would have important implications for R- violating 
operators, since couplings with different family structures would also appear with different 
powers of the family symmetry-breaking parameter. This is consistent with the common 
assumption that a single .R-violating operator dominates. However, this assumption 
would apply in the current quark and lepton basis, and in the mass-eigenstate basis there 
would be several operators corresponding to the original dominant one in the current 
basis. Any given family-symmetry model would make characteristic predictions for the 
pattern of these related operators. Since there are stringent bounds on some of R- violating 
operators, particularly on those involving the first family and on some combinations that 
mix families, an analysis of such sub-leading operators in the mass-eigenstate basis may 
provide the most stringent bounds on the operators related by mass mixing. In addition, 
there could also be further contributions due to operators that are sub-dominant in the 
current basis, with strengths given by powers of the family symmetry-breaking parameter 
that are calculable in any given model. 

The relation between the forms of the mass matrix in the current and the mass eigenstate 
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basis is given by 



K 


= v u L 




M' d 


= v d L 


M ^.(v d R y 


M[ 


= vl 





(1) 

L Ft 

where V udi are the unitary matrices relating the left- and right-handed u, d and t current 
eigenstates to their mass eigenstates. We use the notation C mass = ^f' L M'^/' R for all mass 
terms, so that the Vl are given by diagonalising M'M'\ whilst the Vr are obtained 
by diagonalising M'^M'. Only information on the entries of the Cabibbo-Kobayashi- 
Maskawa (CKM) product matrix 

yCKM = yLiyL (2) 

is provided by experiments to date. 

In general, one can construct models where the quark mixing is either in the up sector, 
or in the down sector, or both. In the class of models studied in this paper, in which 
the mass matrices have small off-diagonal entries generated by spontaneous breaking of 
a family symmetry, one may obtain useful connections between the mixing matrices and 
the elements of the mass matrices in the current basis by perturbative expressions for 
the off-diagonal elements, which are given in the Appendix. From this general analysis, 
it may be seen that in the specific case of the CKM mixing matrix the leading-order 
contribution comes from the ti-quark mass-matrix elements that lie above the diagonal 
in our representation. As a result, we have little experimental input to guide us in 
constructing models for the elements below the diagonal. However, it has has been noted 
for some time that a phenomenologically successful relationship results if one assumes a 
"texture zero" in the (1,1) position and symmetry between the (1,2) and (2,1) matrix 



elements [16]. In this case one finds the relation 




,m d m u m d m u i 

I 1 U 9 . / rn« (h I 2 



V ud \= ( 1 h 2 J COSc 

m s m c y m s m c 



where (f) is the usual CP-violating phase in the CKM matrix. The fact that this relation 
works well is the only phenomenological indication we have for a symmetric structure 
of the mass matrices, and it may just be accidental. Nevertheless, we think it a useful 
starting point for our analysis, so we consider first a simple model capable of yielding 



this form and accommodating the remaining fermion masses and mixing angles [|TT[ . 

The model consists of a single U(l) family symmetry with the same charges for the left- 
and right-handed states, as shown in Table 1, where, e.g., the choice <Zj = (—4, 1, 0) gives 
an acceptable pattern for the mass matrices. Suppressing unknown numerical factors and 
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Qi Ui Di Li Ei H 2 H\ 


17(1) 


d{ di di bi hi — 2a 3 wa 3 



Table 1: Assignments ofU(l) charges. 



phases, which are all expected to be of order unity, with these charge assignments the 
up-quark mass matrix takes the form 



M up = 




The down-quark mass matrix has a similar form, but with a different expansion parameter 
e m y/e. Since the up and down sectors have similar structures, mixing is present in both 
sectors, though it may be larger in the down sector, simply because e > e. For the mass 
matrices of |TT| that we consider initially, one finds the following expressions for the quark 
mixing matrices [) 



/ 1 
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using the second-order perturbation-theory formulae given in the Appendix. 

We now discuss the importance of this mixing for i?-parity violation. The most relevant 
experimental constraints are those on the operators LiQiDi and LiQ 3 D 3 , for which 
A' m < 0.002 from nuclear (3(3 decay ]TJ] and for squark and gluino masses of 200 GeV, 
while A' 133 < 0.001 from bounds [TJ| on Majorana neutrino masses, again assuming masses 
of 200 GeV for the sparticles. Q However, operators related to these operators by mass 
mixing are also strongly constrained by these bounds. Consider first the relations 



{LxQiDx)' = L 1 Q 1 D 1 + eL 1 Q 2 D 1 + 2e 4 L 1 Q 3 £i + •• 
(LiQ 3 D 3 )' = L X QA - zLiQA + e 2 LiQ 3 Di + ... 



2r n n , (3) 



where the notation ()' denotes effective operators in a current-eigenstate basis. We see 
that the operators L1Q2D1 and L1Q3D1 mix with L1Q1D1, so their coefficients are con- 
strained to be the appropriate mixing coefficient (e _1 and (l/(2e 4 ), respectively) times 
the bound on A' m . Similarly, the coefficients of the operators LiQ 3 D 2 and LiQ 3 Di are 
constrained to be less than e _1 and e~ 2 times the bound on A' 133 , respectively. 

We display below, as an example, matrices of upper limits on LiQjDk operators. These 
limits follow from the mixing in this particular model, combined with the experimental 

1 Lepton mixing is discussed in the next section. 

2 The quoted bound is clearly only approximate, as the exact value depends on soft parameters [fT9[ . 
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upper bounds for sfermion masses of 200 GeV |5|, [T7|, [H| |2(], ^T], |2^| . We first look at the 
bounds that arise from mixing with the A' m operator, tabulating the direct experimental 
bounds in cases where they are stronger than those originating from the A' m mixing: 

/ 0.002 0.009 0.04 
0.038(0.009) 0.03 0.3 
V 0.07 0.56 0.001 



T (from 111) 
L ljk < 



In certain entries we have two values, because bounds on eud terms involve mixing in the 
up-quark sector, whereas bounds on v e dd terms involve mixing in the down-quark sector. 
Then we repeat the analysis for mixing with the A' 133 coupling: 



j (from 133) 



0.002 0.04 0.04(0.02) 
0.07 0.03(0.02) 0.02(0.004) 
0.02 0.004 0.001 



and finally we gather all the best limits for the matrix elements in this particular model: 



t best ~ 



I 0.002 0.009 0.04(0.02) 

0.02(0.009) 0.03(0.02) 0.02(0.004) 
V 0.02 0.004 0.001 



TXVV 



Here, the bound on the (2,1) entry arises from constraints on A' 121 from K 
[PT|], in the case that V^^i arises predominantly from the down-quark sector. At this 
stage we have not yet taken into account other bounds, especially bounds on products 
of R- violating couplings that pose even stricter constraints [23J. As an example 
for this Ansatz, the couplings LxQ 2 D\ and L 1 Q 1 Z) 2 appear at such an order in the 
family- symmetry breaking that the strong bound on the product of these couplings from 
contributions to Am K [^3] is not satisfied. We shall return to this and related issues at 
a later stage. 



It should be noted that the model of [jTTJ has mixing in both the up and down sectors. 
Indeed, the (2,3) entry of the down mass matrix is e = 0.23, which is much larger than 
Vg KM . Thus, to obtain viable mass matrices in this example, one needs a suppression of 

■ cos ^5 0. This case, in which the mixing 



f a >rns +a rn^ + 2 , 
v mi, mt V mtmt 



the mixing in | V c b 

between states is much larger than would have been estimated just using the appropriate 
CKM mixing matrix element, serves as a healthy reminder of the potential importance 
of the details of the underlying model for fermion masses when drawing implications for 
i?-violating phenomena. 



3 Exploring Hierarchies of R- Violating Interactions 



We now consider the effect of the U(l) symmetry on the pattern of allowed R- violating 



interactions |10] . We first recall that possible sets of quark and lepton charges leading to 
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correct mass hierarchies are given O] by: 

Case 1: <2j = hi = (—4, 1,0), where ai and hi are the quark and lepton charges respec- 
tively, and 

Case 2: a, = (-4, 1, 0), h = (-f , ±, 0). 

In Case 1, where leptons and quarks have the same charges, one needs an additional 
symmetry in order to eliminate dimension-four nucleon-decay operators. This may be 
done simply by imposing an anomaly- free flavour-independent baryon parity []25|] , under 
which the fields transform as 

Z 3 : (Q, U, D, L, E, H u H 2 ) -> (1, a 2 , a, a 2 , a 2 , a 2 , a) (4) 

This allows only the lepton-number- violating operators, while forbidding baryon-number- 
violating ones f\ 

In Case 2, which is motivated by constaints on HERA- friendly models, the lepton charges 
of the first two generations are half-integers. On might at first think that the residual Z 2 
symmetry of the U(l) forbids the Li 2 QD operators. However, it is straightforward to 
combine this Z 2 with a normal Z^ 1 matter parity, so as to allow these terms while also 
forbidding the UDD terms. This is possible if Z 2 x Z^ is broken to a residual Z 2 by a 
field $ that is odd under both symmetries. In this case, UDD is forbidden, because it 
transforms as (+, — ) under Z 2 x Z^ 1 ■ Similarly, L\ 2 QD transforms as (— , — ) and is also 
forbidden, but it occurs at O (f) through the term ($LQD)/M. 

Let us now pass to the charges of the i?-violating operators. The first thing to notice is 
that the form of the mass matrices only determines the relative charges of the operators, 
not their absolute charges. To see this, note that that the symmetric structure of M up is 
unchanged if we add a family-independent constant to the charges of the U fields. This 
shows that, as we have already mentioned, the charge normalization of our operators is 
undetermined by the mass structure. 

However, anomaly cancellation must be imposed. With the general charge assignment 
given in the first line of Table @, the coefficients of the SU(3) 2 x 17(1), SU{2) 2 x 17(1) 
and Uy{1) x U(l) anomalies are proportional to ^3,2,1, where 

3 3 
^3 = 2 a i + 2 Wl + 2 W2 

11 3 1 

A t = — J^Oi + + -a 3 (w - 2) +4wx + w 2 + 3w 3 (5) 

We demand that these should vanish up to a Green-Schwarz term |[27|| , i.e., A 3 : A 2 : 
Ax = 1 : 1 : 5/3. The effect of this is shown in Table ^ in the first row we have a generic 

3 A flavour-dependent generalisation of this symmetry has been discussed in [p6| . In this case, 
consistent solutions were found containing only a subclass of operators violating lepton number (LLE) 
and baryon- number (U DD). In this way, it was possible to have both lepton and baryon number violation 
without disturbing proton stability. However, we do not pursue such models here. 
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charge assignment where the flavour-independent pieces Wi are to be chosen such that the 
anomaly cancellation conditions are satisfied. Imposing these conditions and reabsorbing 
03 in the definitions of the charges, one obtains the charges shown in the second row of 
Table |2], where a\ = a 



a 3 , and b[ = bi — 63 are of the same form as discussed above, 



i.e., a'i 



-4,1, 



0), 6{ 



-4, 1,0) in Case 1 and b\ 



7 1 

2' 2' 



0) in Case 2. 





Qi 




Li Ei H 2 


H x 


17(1) 


at 




hi h + w 3 -2a 3 


wa 3 


17(1) 




a ■ + iui a ■ — wi 


6- 6- — w\ —w\ 


Wi 



Table 2: Assignments of flavour symmetry charges, before and after imposing anomaly 
cancellation. 

We now discuss the possible hierarchies of i?-violating operators in the two cases. 
Case 1: In this case, the charges of the operators 0^ = LiLjE^ and LiQjD^ are the 
same, and depend only on the values of i,j, k, and not on their order, as given in Table 
0. We note that the constraints on the operators L 1 Q i D 1 from nuclear (3(3 decay and 
on L1L3E3 from bounds on Majorana neutrino masses constrain the choice of the charge 
W\. Since the exact constraint depends on the magnitude of the expansion parameter for 
the i?-violating couplings, we need to consider what the constraints are on this expansion 
parameter. 



ijk 


111 


121 


122 


222 


131 


17(1) 


-12 - wi 


—7 — W\ 


-2- w x 


3 — W\ 




ijk 


133 


333 


223 


233 


123 


17(1) 


-4 - w x 


-Wi 


2-Wx 


1 — W\ 


—3 — W\ 




Table 3: 


Operator 


charges in 


Case 1. 





In the case of the mass matrices, it was suggested in |TTJ that the mixing between Higgs 
fields carrying different U(l) quantum numbers was responsible for filling in the remaining 
elements of the mass matrix. In this case the expansion parameters e and e are as given 
with M 2 , M\ being the mass scales of the heavy Higgs fields H 2 , Hi that mix 



111 



IT 



with the light Higgses responsible for electroweak breaking. The scales of the vacuum 
expectation values < 9 >,< 9 > are bounded from below by (l/Vl927r)M s t r j n9 , the scale 
of the U(l) symmetry breaking. Hence M 2 and Mi are bounded from below by t~ l 9 
and e _1 9, respectively. In the case of R violation, mixing between the operators LLE or 
LQD proceeds through heavy lepton or heavy quark mixing rather than through heavy 
Higgs exchange. If the former are much heavier than the Higgs states, the corresponding 
expansion parameter e' will be much smaller. The limiting case occurs when they have 
string-scale masses, corresponding to 



Mo 



M. 



> 



< 9 > 



string 



M. 



string 



0.02 



(6) 
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Taking this lower limit for the expansion parameter and using the constraint A' m < 0.002 
from nuclear (3(3 decay, we find that | — 12 — w\\ > 2, whilst the constraint A' 133 < 0.001 
from bounds Jl8| on Majorana neutrino masses indicates that | — 4 — W\\ > 2. 

Next, we note that the magnitudes of the couplings in Table |3| are symmetric in the 
three indices ijk. This implies, for example, that at this level the A' 121 and A' 112 couplings 
should have similar magnitudes. This must be made consistent with with the constraint 
(LiQ2Di).(LiQiD2) < 4 • 10~ 9 , which arises from bounds on ArriK [ 23]| . In the present 



context, this constraint indicates that the relevant charge \ —7 — wi\ has to be large, and 
we reach our first HERA-unfriendly conclusion: in this case the e + d — > c interpretation 
of the HERA data would become untenable. 

Thirdly, a related problem is that some couplings to muons would have comparable mag- 
nitudes to those listed in Table §. For example, the magnitude of the A 211 coupling 
would be comparable to that of the A' 121 coupling. However, certain products of cou- 
plings involving electrons and muons have to be extremely suppressed f|. For 200 GeV 
sfermions, 

A231A131 < 2.8 -10- 6 
Kki^2k2 < 3.2 • 10 6 

v lfcl A 2fcl 



A'uiAL < 2-10" 7 



A' A' < 2-KT 7 (7) 



Using the form of the mixing matrices for Case 1: 




we shall see later that these bounds are so severe as to rule out any possible HERA- 
friendly model of this simple type. Note that, in order to obtain correct lepton masses 
within this Ansatz, a factor of ~ 3 is needed in the (22) element of the mass matrix, and 
this factor also enters in the mixings. 

Other strong constraints on products of couplings are the following: 

AijiAij2 ^ 2.8 ■ 10 ^ 

a: 13 a: 31 < 3.2 -io- 7 

a: 12 a: 21 < 4-io- 9 , (8) 

which are particularly stringent in the model under consideration. Using these bounds, 
one finds 

Ka < 6 • IO" 4 

K12 < 6 • 10~ 5 , (9) 
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together with corresponding bounds for permutations of the indices. 

If the reported apparent excess of HERA events at large Q 2 were due to production of a 
single squark by an i?-violating coupling, one would need 

A' m , 13 i~ 0.04/ VB 
A' 132 ^ 0.3/VB 

where B is the branching ratio of the decay q — > e + q. We see immediately from @ that 
the first two possibilities cannot be realised in this model. Moreover, we see from ([3D 
that A' 132 ~ A' 133 /e < 0.004, and infer the third possibility cannot be realised either. 

What happens to the remaining couplings? From the mixing discussed above, we have: 

(i) A' m = A' 112 /e < 3 ■ 10~ 4 , which is a stronger bound than the one from neutrinoless (3f3 
decay, Q 

(ii) A' 222 = A 221 /e<3-10- 4 , 

(iii) A 223 = A 213 /e = A' 312 /e < 3 ■ 10~ 4 or A' 223 = A' 213 /e = A' 312 /e < 0.0013. 

The constraints here are very strict because the experimental bound on A 312 is more 
severe than that on A 213 : since we require these two terms to have the same charge, we 
must take the stricter limit. We also have 

(iv) A' 233 = A 223 /e < 0.0013(0.006), for expansion parameters e and e respectively, 

(v) A' 333 = A 233 /e< 0.006(0.025). 

In each of these cases, R violation may be manifest in hadron-hadron colliders. For 
sfermion masses of 100 GeV and A > 10~ 6 , the lightest supersymmetric particle is ex- 
pected to decay inside the accelerator. The above constraints allow couplings that are 
significantly larger than this lower bound. Through a suitable choice of w\, the couplings 
that are more severely constrained can be made small, while some others can be of im- 
portance for collider physics, though none can be very large in this type of model with 
symmetric mass matrices. Hence, single-squark production via an i?-violating coupling is 
suppressed, and the best signal would be squark-pair production followed by _R-violating 
decay. 

Case 2: 

In this case, the charges of the operators depend on the flavour-symmetry charge of the 
singlet field $ that we have introduced. This does not affect the relative magnitudes 
of the i?-violating couplings, since $ appears in all terms. However, this charge and 
the vacuum expectation value of $ do provide a possible source of suppression for the 
R- violating couplings. We take as an indicative value a$ = 1/2: the corresponding 
subclasses of LLEQ and LQD<& operators with integer flavour charge appear in Tables |] 
and |5|. 

In this model, the lepton mass matrix takes the form 

4 We use the down-quark mixing parameter, which is larger, since the operator we compare with differs 
only in the index of D. 
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ijk(LLE) 


121 122 133 233 


17(1) 


—6 — w\ — 2 — wi — 3 — Wi 1 — Wi 



Table 4: Integer LLE charges for Case 2. 



iHLQD) 


111 


121 


122 131 123 


133 


U{1) 


-11 - w 1 


—6 — wi 


— 1 — wi —7 — Wi —2 — wi 


—3 — wi 


ijk(LQD) 


211 


221 


222 231 223 


233 


17(1) 


—7 — W\ 


-2-w x 


3 — W\ — 3 — u>i 2 — wi 


1 — 10i 



Table 5: Integer LLD charges for Case 2: these charges remain the same if j and k are 
interchanged. 






which is independent of the value chosen for w±. What are the predictions for the strength 
of the i?-violating couplings in this model? As in the previous model, we have: 

Kis < 6 • io- 4 

Ku < 6 ■ io- 5 (10) 

and so again there is no possibility to explain the HERA events, essentially for the same 
reasons as in Case 1. From the charges of Table ||, we find that A^n = A' 113 , A' 133 = A' 213 
and A' 123 = A 212 - Since A' 113 = A' 123 /e < 3 ■ 10~ 4 , we obtain the slightly stronger bound 
Kn = K31 — 3 ■ 10~ 4 . For the remaining couplings we have the following bounds: 

(i) A' m = A' 112 /e < 3 • IO" 4 , 

(ii) A' 122 = A' 121 /e<3-10- 4 , 

(iii) A 222 = A 221 /e<3-10" 4 , 

(iv) A' 223 = A 213 /e < 0.0013, 

(v) A' 233 = A 223 /e< 0.006. 

The difference from the previous solution is that the couplings L^QjDk are absent, and 
thus constraints from them are evaded. However, the model remains restrictive, as all 
the quarks of the same generation have the same charge. Therefore, the strict bounds on 
products of operators still constrain strongly individual couplings. Nevertheless, we see 
from the limits above that several possibilities exist for .R-violating squark decays within 
hadron-hadron collider detectors. 

We see therefore that there are four problems that do not allow an explanation of the 
HERA events within the framework of these models. First, the quarks and leptons of 
the same generation have the same charges, so the LiLjE k and LiQjD k couplings are 
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subject to the same bounds. Secondly, the choice of symmetric mass matrices makes 
the last two equations of (|8]) difficult to satisfy, because Qi and Di have the same charge 
and hence each factor involved has the same suppression, so that one cannot arrange to 
satisfy the inequality while keeping one coupling large. Thirdly, the model has large 
mixing in the (1,2) down-quark sector, making the last of eqs (H) difficult to satisfy. 
Finally, the large (1,2) mixing in the charged lepton sector may not be reconciled with 
bounds on products of couplings that involve electrons and muons. 

Thus we see that the combination of the various i?-violating bounds with simple family 
symmetries produces strong constraints on a variety of i?-violating couplings. For the 
case of the family symmetry leading to the symmetric mass matrix (Q) these constraints 
imply that R violation does not give rise to anomalous events at HERA at a significant 
rate. 



4 HERA- Friendly Textures of R- Violating Couplings 

In this section we explore modifications of the simple U(l) family structure, which may 
be able to accommodate an .R-violating interpretation of the apparent excess of events 
at HERA. As we have stressed, a major problem in building a model to accommodate 
the HERA events lies in the need to satisfy the bounds (|8|) while keeping large one of 
the individual couplings involved in these products. This leads us to consider models 
with the (1,2) mixing entirely in the up-quark sector and to deviate from the symmetric 
mass-matrix structure. 



4.1 Asymmetric Flavour Textures 

Once one gives up on the symmetric form, the pattern of masses is insufficient to constrain 
the U(l) charge structure, so there are many new possibilities. Here we present just one 
viable choice to illustrate the options, but we certainly do not claim any uniqueness. We 
start with the charge assignment (-4,1,0) for the quark doublets of the model discussed 
above, and modify the up- and down-quark singlet charges to achieve the desired struc- 
ture. In order to reduce the arbitrariness, we also choose to generate both the up- and 
down-quark mass matrices with the same expansion parameter, as would be the case if 
the non-renormalisable terms QiDiH20/M arise through heavy-quark mixing. 

With this Ansatz, a suitable choice for the up-quark singlet charges is (-5,1,0), which 
gives 

/ e 9 4e 3 £ 4 \ 

M up = 4e 4 e 2 e (11) 
U 5 e 1 ) 

A value of e ~ 1/20 gives an acceptable charm mass. Note we have been forced to assume 
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an enhancement factor of 4 in the (1,2) and (2,1) elements in order to accommodate the 



mixing needed to generate V^ 2 ~ = 4e ~ 0.2. The mass eigenvalues are 1, e and 
16e 5 , from which we see that the factor of 16 which appears from the coefficients in the 



off-diagonal entries compared to the solution of [11] is compensated by the additional 
power in the expansion parameter. 

The choice of down-quark charges is dictated by the requirement that we keep the (1,2) 
mixing small. A suitable choice for the charges of the singlet down quarks is (7,-3,1), 
which gives the structure 



j^down 



( e 3 e 7 e 3 

e 8 e 2 e 2 I (12) 



\ e 7 e 3 



The eigenvalues scale as e, e 2 ,e 3 , and V% ~ V^ zi ~ e 2 . Moreover, VP ~ e = 0.05, 
so we do not require the cancellations that were needed in Jl l| (remember that in this 
case V^ 3 32 ~e = 0.23). Note that this choice has the advantage of reducing the bottom 
mass through an e factor, putting us in the small-tan/3 regime. This phenomenological 
choice of charges does not yet ensure anomaly cancellation, but at a later stage, when 
we also have a good phenomenological choice for the lepton mass matrix, we will discuss 
what flavour-independent charges have to be added in order to cancel anomalies. These 
additional charges will not modify the hierarchy of couplings, nor the relative magnitudes 
of the i?-violating couplings of a given type. 

The key point of this model is its large U(l) charge difference between the relevant LQD 
couplings: 

a L 1 Q 2 D 1 - a L 1 Q 1 D 2 —* 15 (13) 
O^OsSi - ^iQiA, 10 ( 14 ) 

which leads to large relative suppressions of these operators, though mixing in the (1,2) 
and (1,3) down sectors will close this gap. Consider first the operators appearing in (pf). 
The mixing of the left-handed down quarks is given by the form of 21 in the Appendix, 
and the second-order term dominates with rnf 3 mf 2 / {M3M2) ~e 3 = 1.3- 10~ 4 . The mixing 
of D\,r>2 is given by m^/il^ ~ ^31^23/(^3^2) ~ e 6 - Taking the same expansion 
parameter as for the masses, consistent with @, the net suppression is e 3 e 6 ~ 2 ■ 10~ 12 . 
This is more than sufficient to satisfy the bound on L1Q1D2 while allowing the L1Q2D1 
to have a coefficient large enough to give the HERA events. Similarly, one may check 
that it is possible to for the L1Q3-D1 operator to be relevant for HERA, without inducing 
an L1Q1D3 coupling with an unacceptable value. 



4.2 Lepton Flavour Violation 

We now turn to the assignment of lepton charges in such a model. Given the bounds of 
([7D, if any coupling involving an electron is large, the corresponding coupling involving 
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muons should be small. The simplest solution is to choose charge assignments so that 
the U(l) charge of the (1,2) entry of the lepton mass matrix is half-integer. In this case, 
a residual Z 2 symmetry forbids the (1,2) mass term. The choice of charges is restricted 
by the anomaly cancellation conditions. These are: 

A 3 = (a Ql +a Q2 +a Qa ) + ^(a 0l + a 0i + a A + ^(a Dl + a B2 + a D A 

3 11 
A 2 = -(a Ql + a Q2 + a Qs ) + -(a Ll + a L2 + a Ls ) + -(a Hl + a H2 ) 

1 4 1 

M = ~(a Ql + a Q2 + a Qs ) + -(a Dl + a 02 + a 03 ) + -(a Dl + a D2 + a Ds ) 

+^i a L 1 + a L2 +a L:j ) + (a Sl + a Si + a &i ) + -(a Hl + a H2 ) (15) 



3 



where by we denote the charge of particle F in the i generation. We see from A 
and A 2 that, for integer quark charges, a natural solution of the conditions has the sums 
of a#. and integers. Thus, we need two of the ag. and to be half-integers. We 
see for example, that the choice a Ll = 9/2, a L2 = — 1, Ol 3 = —1/2,0^ = —l/2,a^ 2 = 
— 1, a^ 3 = —1/2 generates a viable mass hierarchy. Here we have chosen the (3,3) charge 
to be the same as that of the down-quark, in order to give b — r unification. Clearly this 
pattern of charges is not the only viable choice, but it indicates how things may work. 
The corresponding lepton mass matrix in this solution is: 



M 



( e 



\ 



e 4 





e 



(16) 







and there is (1,3) mixing, but no (1,2) or (2,3) mixing. The eigenvalues of this mass 
matrix are e, e 2 , e 4 , consistent with the measured values for e = 0.05. 





Qi 


Ui 


A 


Li 


Ei 


H 2 Hi 


U(l) 


a! t 


a 'u, 











U(l) 


a'i 


a 'u, + w i 


—1+ a' D . - wi 






— W\ 1 + W\ 



Table 6: Flavour charges in models with asymmetric mass matrices. 



We denote by a' Qi = (-4, 1, 0), d % = (-5, 1,0), a' Di = (7, -3, 1 



' Li = (9/2,-1,-1/2) 

and a'g = (—1/2, —1, —1/2) the choices of generation-dependent charges in this model. 
The top line of Table |6], which includes these and the corresponding charges for the 
Higgs multiplets, is not anomaly-free. It is easy now to satisfy the anomaly- matching 
conditions, allowing for additional family- independent components of the U(l) charge, 
which do not affect the mass matrix structure. An anomaly-free solution is obtained by 
adding charges as indicated in the second line of Table |6|, where the variable wi is an 
integer. 



13 



4.3 Nucleon Stability 



We now demonstrate that nucleon decay graphs due to combinations of LQD and U DD 
interactions may be eliminated in this HERA-friendly example by imposing an anomaly- 
free discrete gauge symmetry. In the specific model discussed above, where the first- 
and third-generation leptons have half-integer charges under the flavour symmetry, we 
can again combine the residual Z2 symmetry of the U(X) with a normal Z^ 1 matter 



parity, where Z 2 x Z^ 1 is broken to a diagonal Z 2 by a field $ that is odd under both 
symmetries. Then the couplings UDD and L2QD transform as (+, — ) under the sym- 
metry and are forbidden. Renormalisable L\^QD couplings which transform as (— , — ) 
are also forbidden, but effective couplings of this type may occur at O [$^\ through the 
term ($LQD)/M. 

Such underlying Z2 symmetries can be illustrated in the context of a GUT group, if de- 



sired. Consider, for example, the Pati-Salam gauge group SU (4) x SU (2) L x SU (2) R [^8 



In models based on this group, the fermionic fields belong to either the 4 or the 4 repre- 
sentations of SU(4), and no trilinear R- violating term is invariant under the symmetry. 
However, invariants can be constructed by introducing an adjoint field £ [p9[| . If $ has 
a half-integer charge, the fact that it is in the adjoint of SU(4) means that all baryon- 
number-violating operators are forbidden in any order, whilst the terms Li 3 QD§ have 
integer charge and are therefore allowed. Moreover, no effective terms Li 3 EH 2 Q, which 
could cause problems with the lepton mass hierarchies, are allowed, as they are not 
invariant under the extended gauge group. 



4.4 Hierarchy of R- Violating Interactions in a HERA-Friendly 
Model 

We now consider the effect of the U(l) symmetry on the pattern of allowed R- violating 
interactions in the models that were motivated by the HERA events. The charges of the 
operators depend on the half-integer charge of the field $ under the flavour symmetry. 
This does not affect the relative magnitudes of the R- violating couplings, since $ appears 
in all terms. However, this charge and the vacuum expectation value of $ do provide a 
possible source of suppression for the R- violating couplings. The corresponding subclasses 
of LLE§ and LQDQ operators with integer flavour charge, before introducing mixing 
effects, appear in Tables [7| and We have used the charges of Table imposing anomaly 
cancellation. Here we have taken <3> to have U(l) charge 1/2, but its actual value can be 
re-absorbed in the definition of W\. 

Let us first consider the possibility that the possible excess HERA events are due to the 
L1Q3D2 operator with A' > which will be the dominant operator if w\ = 0. 

The relative suppression of the L1Q3D3 operator is e 4 . Of course, mixing effects in the 

5 A study of the string origin of non-renormalisable operators in this model has been presented in j30[. 
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(2,3) sector of the down-quark mass matrix re-introduce an L1Q3-D3 operator at order 
e 2 , via the mixing of right-handed down quarks. In the type of models with a small 
Vj£ mixing therefore, this solution may in principle be accommodated. However, for 
our specific choice of charges, we see that the operator LiQ s D 2 has the same charge as 
L1Q1-D3, which is bound by charged- current universality [^(J to be < 0.04 for a squark 
mass of 200 GeV. Hence the L1Q3D2 interpretation of the HERA data is not realisable 
in this specific model. However, this is rather accidental for this particular example, and 
need not be the case in general. 



ijk(LL£$) 


122 131 133 231 


17(1) 


1 — W\ 2 — w\ 2 — w\ —4 — w\ 



Table 7: Integer LLE& charges, ignoring mass mixing. 



ijk(LQD<5>) 


111 


112 


113 


121 


122 


123 


U(l) 


6 — Wi 


-4 - w 1 




11 -Wt 


1 — Wt 


5 — ^i 


ijk(LQL>$) 


131 


132 


133 


311 


312 


313 


17(1) 


10 -w x 


-Wi 


4- wt 


1 — Wt 


-9 -wt 


—5 — wt 


ijk(LQ7J$) 


321 


322 


323 


331 


332 


333 


U(l) 


6 — Wi 


-4 - wt 


-Wt 


5 — Wt 


—5 — wi 


—1 — Wt 



Table 8: Integer LQD& charges, ignoring mass mixing. 

We now look at the possibilities that the HERA events arise from the LiQ2,3-Di operators. 
We see from Table |6] that, in the absence of mixing, the relative suppressions of the 
LtQtDt and L1Q3-D3 operators would have been enough to make these cases viable. 
When mixing effects are included, the possible effects of unknown phases should be taken 
into account when comparing with bounds. In the specific case that the HERA events 
are due to an L1Q2-D1 coupling, we have no problem with the LiQ 3 D 3 operator, but 
there is a potential difficulty with (3(3 decay, due to mixing with the LtQ\Dt operator: 



'200 GeV s 



< 4 • 10" 



m- n 



1 TeV 



1/2 



;i?) 



where nig is the gluino mass. Given (a) that the V mixing arises from the up 
sector in our framework 0, and (b) that the bounds from the Tevatron indicate that 
the branching ratio of cl to fermions can not be close to unity in the context of this 
interpretation, implying that A' 121 has to be larger than 0.04, we see that this solution is 
not naturally accommodated. It might be possible if m>u L is significantly larger than m>c L , 

6 Even in the case that the V^^^ 1 mixing arises from the down-quark sector, squark mass universality 
violation is required in order to evade bounds from K — ■» irvv. 
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but this requires a violation of squark-mass universality that is potentially dangerous for 
flavour- changing neutral interactions. 

On the other hand, if the HERA events are due to a L\Q 3 D\ coupling, there is no problem 
with the L1Q1D1 operator, but a problem could in principle appear with the LiQ 3 D 3 



coupling that is bounded from limits on neutrino Majorana masses However, the 
relevant (3,1) mixing term is small, indicating that in this case the unknown coefficients 
may be such that the bounds are easily accommodated. What about the other couplings? 
For LiQ 3 Di « 0.04, the model predicts that L 1 Q 2 D 1 m L 3 Q X B 2 « 0.002, while all other 
couplings are very suppressed. Indeed, looking at the charges, we see that the next 
larger couplings are suppressed by e 4 as compared to LiQ 3 D\. Mixing effects are also 
suppressed, except for the operators L\Q\D\ <m L\QiD\{4d) ~ 0.0004 and L 3 Q 2 D 2 ~ 
L 3 QiD 2 (4e) ~ 0.0004, which are within the allowed range. 

Finally, note that we do not have any mixing between L\QjD\. and L 2 Q^D^ couplings 
(the later are forbidden by the symmetry), so the dangerous product combinations that 
violate lepton flavour are also absent. 

In the light of the above discussion, we conclude that, of the valence quark produc- 
tion mechanisms via LiQ 2 D 1 and LiQ 3 Di couplings, the second possibility seems to be 
favoured. It should be possible to make a model with a coupling L\Q 3 D 2 sufficiently 
large to explain the HERA data, although we have not displayed one here. 



5 Baryon Decay via Dimension-Five Operators 

We saw earlier on that the experimental absence of baryon decay imposed important 
constraints on possible models, which are most easily evaded by imposing a baryon parity 
symmetry that forbids the dangerous UDD couplings. However, this is not the end of 
the story, since models may also contain dimension-five operators that would generate 
proton decay at an unacceptable level. The most dangerous among these operators are 
the operators [QQQL}p and [QQQHijp, the latter in the presence of LQD couplings. [] 
These operators can lead to fast proton decay via loop diagrams. In the case of [QQQL] F 
operators that involve the two lightest generations, the constraint on the coupling rj of 
any such operator is rj < 10~ 7 . This bound has some flexibility, since the magnitude of 
the loop diagrams depends on details of the sparticle spectrum, but this possibility is 
not crucial for the subsequent discussion of models. In the case of [QQQHi]p operators 
with couplings rj', fast proton decay may occur if they are present simultaneously with 
LQD operators with generic coefficients A'. The product of the corresponding couplings 
is constrained: rj'X' < 10~ 10 . Since an /^-violating interpretation of the HERA events 
requires either LiQ 2 Di or LiQ 3 Di « 0.04, it is clear that we have to worry about the 
[QQQHi}p operator as well. 

7 The lepton-number-violating operators [QUEH±]f, [QUL*]d and [QUL*]d are dangerous in the 
presence of UDD ones. 
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We now analyse the dimension-five operator charges in the different cases discussed in 
previous sections, to see whether they are large enough for the suppression by powers 
of small quantities to be sufficient. How small the terms actually are depends on the 
expansion parameter, as we have already discussed in a previous section. 

The QQQH\ operators are easily dealt with, even though the baryon stability require- 
ments seem to be more severe for them. The reason is that these operators transform 
as (+, — ) under the Z 2 x Z 2 l , and are thus forbidden. What about the QQQL opera- 
tors? The QQQLi s operators are not present, because they transform as (— ,+) under 
Z 2 x Z^ 1 . However, the operators QQQL 2 are allowed. These are dangerous, because 
proton decay may occur via the modes p — > ^1,2,371" + and p — > v\^ 2 ^K + . Let us look at the 
flavour charges of these operators. We recall that colour antisymmetrisation implies that 
all the quark flavour indices cannot be identical. The operators that are not suppressed 
enough by quark mixing parameters have the following charges in the model that could 
explain the HERA events: 

a QiQiQ2L 2 = -9 
a QiQ 2 Q2L 2 = -4 
a QiQiQ3L 2 = -10 

a QiQ2Q:iL 2 = -5 (18) 

where for the lepton charge we used the anomaly-free choice of Table ^. 

We infer that we do not need any further underlying symmetry in order to suppress these 
couplings adequately. However, even in models where this suppression does not occur, 
there could be some GUT symmetry that forbids the offending QQQL operators 0. This 
would be an interesting constraint on GUT model-building, but should not be taken as 
a serious obstacle to constructing HERA-friendly models. 



6 Concluding Comments 

We have discussed the implications of a single U(l) abelian flavour symmetry for the 
possible hierarchies of i?-violating couplings. The relations between the Standard-Model 
Yukawa couplings and i?-violating couplings depend on the choice of model charges, so 
the observed hierarchies of quark and lepton masses do not lead to a unique specifica- 
tion of the dominant /^-violating couplings. However, we have identified certain general 
features of such a framework, highlighting the importance of mass mixing between cur- 
rent eigenstates. We have identified various interesting possibilities for hadron-hadron 
collider phenomenology that are consistent with this mixing and the available experimen- 
tal constraints. Within this general approach, we have searched specifically for simple 

8 Moreover, in string-derived GUT models, string selection rules may lead to the vanishing of operators 
that are invariant under field-theory symmetries. In such models, it is possible to construct realistic 
fermion mass matrices while having maximal proton stability |fsi| . 
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consistent models that lead to the favoured i?-violating scenarios for the explaining the 
possible excess in the HERA data. 

Our results may be summarised as follows: 

• Flavour symmetries lead us to expect a hierarchy in the R- violating couplings, analogous 
to that observed for the known fermion masses. These hierarchies can be consistent with 
a squark- production interpretation of the HERA data (if required), as well as with the 
various other experimental constraints on the couplings. 

• The simplest charge assignments lead to unified, and thus more predictive, forms for 
the mass matrices. For the case of equal charges for up and down quarks and leptons 
of a given generation, the symmetry together with bounds from products of i?-violating 
couplings implies that there should be no significant anomalous events at HERA coming 
from such couplings. If we wish to accommodate such anomalous events, we are forced 
to depart from this picture. Schemes with asymmetric charges and different assignments 
for up quarks, down quarks and leptons give rise to larger splittings between different 
operators. 

• Some of the charge assignments considered forbid large coefficients of dimension-five 
operators that are potentially dangerous for baryon stability. In schemes where this is 
not true, such terms would need to be forbidden by further GUT symmetries. 

One can consider relaxing various of our conditions, for example by introducing a higher 
level of asymmetry in the mass matrices, invoking multiple U(l) flavour symmetries, 
etc., and in such models the predictions can be further altered. Moreover, additional 
zero couplings may be expected when one goes to a specific GUT/string construction. 
However, it is interesting that it is possible to construct phenomenological models with a 
single U(l) flavour symmetry that are compatible with attempts to explain the reported 
excess of HERA data by i?-violating squark production, albeit at a price. In order to 
constrain the possible schemes, and perhaps rule some out, more experimental data are 
required. 



Appendix 



Using second-order perturbation theory, it is easy to derive the mixing elements for a 
generic mass matrix 



52[ |33|, where m stands for the off-diagonal contributions and M 
for the diagonal part. The left-handed mixing is given by |33| 

(m ik M k + mliMJimuMj + m* k M k ) 



VA 



Mf 



Mj 



+ 



m ik m jk 



Mf, 



(Mf — Mj)(M% — Mj) (Mf 

and V$ is given by a corresponding expression, substituting the mass matrix by its 
hermitian conjugate. In the case that the ratio of rriij to m*j is considerably larger than 
the ratio M^/Mj, the mixing elements are given by: 



U, 
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' Mo 



m 13 m 32 
M^Mo 
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+ 



m 32 m 13 
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M 2 M 3 
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12 Mt 23 
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M 2 M 3 
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+ 



(19) 



In the above expressions, terms in brackets mark contributions which involve mass entries 
below the diagonal. These terms, as well as the ones in parentheses, can in most cases be 
neglected. However, they may become relevant if the mass matrices have texture zeroes. 

In the case of the full Vqkm matrix, one has 
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(20) 



These formulae are used in the text in conjunction with specific parametric forms for the 
off-diagonal terms in the up- and down-quark mass matrices. 
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